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(-H I We present a first status report on a broad and systematic study 

of possible manifestations of QCD-instantons at HERA. Considerable 
motivation comes from the close analogy between instanton-induced 
^^ . B + L violation in electroweak processes and effects of QCD-instantons 

5h ' in deep inelastic scattering. We concentrate on the high multiplicity fi- 

^ nal state structure, reminiscent of an isotropically decaying "fireball" . 

A set of experimental isolation criteria is proposed. They serve to 
further enhance the striking event signature without significantly sup- 
pressing the expected rates. 
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1 Introduction 

The basic significance and possible importance of QCD-instanton effects in 
deep inelastic scattering for decreasing Bjorken variable x-qj and high photon 
virtuality Q^ has recently been emphasized 0]. 

First of all, it has been argued pi M that the calculable small scale 
instanton dynamics [^ may be factorized from the large distance effects, 
allowing for a semiquantitative estimate of instanton-induced contributions 
to the structure functions ||l| and the hadronic final state Q. This is rooted 
in the fact that Xj \fQ^ plays the role of a dynamical infrared cut-off for the 
instanton size p. 

Secondly, QCD-instanton effects for decreasing xgj are largely anal- 
ogous to the manifestation of electroweak instantons for increasing energies 
[y] . The anomalous B + L violation due to electroweak instantons is paral- 
leled by a chirality violation induced by QCD-instantons |^. The spectacular 
possibility of a strong increase of ranlti-W^ / Z^ production in the multi-TeV 
regime |^, ^] due to electroweak instantons corresponds to a strong enhance- 
ment of multi-gluon production at small xbj due to QCD-instantons. Strik- 
ing consequences include a high multiplicity final state structure, reminiscent 
of a decaying "fireball" . 

Whereas a promising search for anomalous electroweak events is only 
possible in the far future, presumably at a post-LHC collider |8| or at cosmic 
ray facilities Q, the search for anomalous events induced by QCD-instantons 
can start right now, in deep inelastic e^p scattering at HERA. 

The present paper represents a first status report on a broad and sys- 
tematic study of possible manifestations of QCD-instantons which could be 
searched for at HERA [Q]. While a theoretical derivation of crucial quanti- 
ties characterizing the instanton-induced final state (1,2... particle inclusive 
rates, average multiplicity, average transverse momenta, energy flow, etc.) 
will be deferred to Ref. ^, we mainly report here on the phenomenological 
aspects of our results. 

The organization of this contribution is as follows. We start off 
in Section 2 by expanding on the close analogy between instanton-induced 
B + L violation in electroweak processes and effects of QCD-instantons in 
deep inelastic scattering. Section 3 contains a summary of the results of Bal- 
itsky and Braun ^ concerning the QCD-instanton induced contribution to 
the gluon and quark structure functions. We also emphasize the approxima- 
tions and limitations inherent in this calculation. Furthermore, we present 
the instanton contribution to the nucleon structure function F2{x-Q-pQ'^) 



obtained by convoluting the results of Ref. [^ with phenoinenological distri- 
butions of quarks and gluons in the nucleon. In Section 4 we report on our 
ongoing investigation 0] of the QCD-instanton induced hadronic final state 
in deep inelastic scattering. The main emphasis rests on the characteristic 
event topology along with a discussion of experimental isolation criteria, 
serving to further enhance the striking event signature without significantly 
suppressing the rates. A search strategy for instanton- induced events is for- 
mulated. In Section 5, we present a summary and an outlook on related 
aspects and open problems under study. 



2 The QFD - QCD Connection 

The Standard Model of electroweak (QFD) and strong (QCD) interactions 
is remarkably successful. In particular, its perturbative formulation ("Feyn- 
man diagrammatics" ) appears to be theoretically consistent and agrees with 
precision experiments (where applicable, i.e. for small coupling constant). 

Nevertheless, even for small couplings, there exist physical processes 
which cannot be described by conventional perturbation theory, notably, 
phenomena associated with quantum tunnelling. 

In non-Abelian gauge theories such as QFD and QCD the vacuum 
actually has a complicated structure, even on the classical level |jl^: The po- 
tential energy is periodic with respect to the so-called Chern-Simons number 
(c.f. Fig.|lD, 

A^cs [^] = yI^ / d'x e,,u [Atd.Al - ^-ea,cAtA]A%) , (1) 

which is the (topological) winding number of the (non-Abelian) gauge field 
[A) under consideration. Pure gauge fields corresponding to the degenerate 
minima of the potential energy (perturbative vacua) have integer values of 
the Chern-Simons number A'^cs- Moreover, pure gauge fields differing by 
AAcs = n = integer are topologically inequivalent. They are related to 
each other by a topologically non-trivial static gauge transformation with 
winding number n. This means that they are separated by an energy barrier, 
as shown schematically in Fig. ||. 

In the electroweak theory a mass scale, v ~ 246 GeV, is introduced 
via spontaneous symmetry breaking and, correspondingly, there is a definite 




Figure 1: Schematic illustration of the static potential energy of the gauge (and 
Higgs) field vs. the Chern-Sinions number 



minimum barrier height associated with the VF-mass, of order ||Tl| 
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As is well known, this minimum barrier energy is associated with a cer- 
tain static, unstable solution to the classical field equations, the so-called 
'sphaleron' |11|. It may be viewed [12| as an intermediate, coherent field 
configuration consisting of a large number, 0{l/aw)-, of W/Z (and Higgs) 
bosons, confined to a characteristic volume of order m^ . 

In QCD, being an "unbroken" gauge theory, the minimum barrier 
height depends on the considered process and its kinematics. It turns out, 
that the notion of the minimum barrier may be qualitatively transcribed 
from QFD, provided that there is a (process dependent) hard scale Q' avail- 
able. Specifically, in deep inelastic e^p scattering. 



Q'^ oc g2 



-(mass) of the virtual photon. 



(3) 



Besides aw -^ ««, one is then led to substitute the characteristic linear 
dimension ||, |, 0, |l3| 

— ' ^ • (4) 

mw as{Q')Q' 

The minimum barrier energy ( from Eqs. (§,§)) 



<a™r ~ Q'^ 



(5) 



is now associated with a sphaleron-like intermediate, coherent field config- 
uration consisting of a large number, 0{l/as), of gluons in a characteristic 
volume {asQ')~^- 

Transitions between minima of the effective potential in Fig. |l| lead 
to a violation of fermionic quantum numbers in the Standard Model [^, ^ . 
In particular, baryon [B) and lepton (L) number conservation is violated 
due to non-perturbative electroweak gauge fields {W) associated with the 
SU(2) flavour gauge group, according to the selection rule 

ALe,^,. = ^A5 = -AiVcs[VF]. (6) 

In analogy, non-perturbative gluon fields {G) associated with the SU(3) 
colour gauge group induce a violation of chirality conservation for (massless) 
quarks, 

AQ5^,^,,... =2AiVcs[G]. (7) 

For parton-parton center of mass (cm.) energie^ vV < Mbarrier, 
such processes are classically forbidden and only occur via quantum tun- 
nelling under the barrier in Fig. |l|. In this case, the sphaleron-like inter- 
mediate state corresponding to the energies (y), (|5|) can only be reached 
virtually. Correspondingly, the respective cross sections are exponentially 
suppressed in the coupling [^, 

o-tunneiiing ~ exp(-47r/a) ; a = awiOLs- (8) 

Let us recall the origin of this tunnelling suppression factor in some- 
what more detail. The amplitude of anomalous fermion number violating 
processes can be obtained by expanding the path integral for the correspond- 
ing Green's functions about instantons [15, ^, 16|, which are classical solu- 



tions of the Yang-Mills(-Higgs) equations in 4-dimensional Euclidean space 
with finite action. The instanton (/) (anti-instanton (/)) may be viewed as a 
most probable tunnelling solution, interpolating in Euclidean time between 
the gauge (and Higgs) field configurations of two neighbouring vacua with 
AA'^CS = +1(~1)- It passes the above-mentioned sphaleron-like intermedi- 
ate state in-between. 



^ Henceforth, we shall denote by y/W the total cm. energy of the non-perturbative 
subprocess, while -^s refers to the total cm. energy of the physical process under 
consideration. 



In QCD, for example, the instanton is explicitly given by]^ 

i Ui[a^{x - xi) - {x- xi)^]Ui 



G^^\x-xi-Ui,pi) 



{x - xiY{{x - xiY + pf) 



pI 



(9) 



depending on a set of collective coordinates {x/,/9/, C//}, like center, xi, 
size, Pi, and orientation in group space, Uj. Since the action is independent 
(QCD) or only slightly depends (QFD) on these collective coordinates, they 
are to be integrated over. For the simplest exclusive anomalous processes 
one obtains (for QFD, see Refs. [^, ^, for QCD, see Refs. [§, H, ||]) in this 
way to exponential accuracy: 



A{B + L) = -2nfamiiy = -6 : 

T [qq -^ 7q3l) oc / dpj ... exp 
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(11) 

where, in our normalization, the action of the QCD-instanton is given by 



^QCD 



^QFD 



1 + 



1, whereas the action of the QFD-instanton reads Sj 
(l/2)/9j?Ti,^. In the various QCD-instanton induced subprocesses of Eq. ([ll|), 
the star indicates that the corresponding parton carries a virtuality Q' = 
\/—q^ > 0. We observe that, for exclusive anomalous processes, the instan- 
ton size is effectively cut off at pcut ~ v~^ in the electroweak theory ||5|, |^] 
and at pcut ~ Q'^^ in QCD |^ |2|, |l[ (for early discussions, see Ref. p!7|]). 

Most interesting from a theoretical point of view and in the light of 
present and future collider experiments is the case \/s' > Mbarrier, where a 
transition over the barrier is classically allowed, i.e. energetically possible. 
Unfortunately, the crucial dynamical question is still unsettled |^ whether 
the transition from a state with a few initial partons to the very different 
sphaleron-like multi-parton coherent state can proceed without extra sup- 
pression. Only then could such anomalous processes acquire an observable 



^ We use here the notations a"" — {—ia, 1), c 
matrices). Furthermore we abbreviate x — x^a^ 



— {+i<J, 1) (f are the standard Pauli 



, X — X uG ij 



etc. 



cross section and the final state would consist of a large number, 0{a^^,^^), 

of W/Z^s (gluons) in addition to the few fermions required by the anomaly. 

The possibility that this intriguing scenario might be realized in 

nature was first observed in the context of the electroweak theory in Ref . [0 . 



It was found |18| that, to exponential accuracy, the total cross section for 
anomalous B -\- L violation, in the high energy and weak coupling limit. 
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aw 



fixed, 
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can be written in the following scaling form 



(I) tot 
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tx exp 






(12) 



(13) 



where Mq = vGvrmiy/avy is of the order of the minimum barrier height Eq. 
(0). The so-called "holy-grail" function F in the exponent is known only in 



a low energy expansion whose first few terms are given by [19| 

FQFD(e) = 1 - ^ 6^/3 + ^e' + 0{e'/^[l + In 6]), 
8 lb 



(14) 



where e = Vs'/Mq. The effective instanton size entering the running cou- 
pling in (13) scales like m^ , 



1 



mw 



:2/3 



(15) 



Note, that the first term in the series expression for the holy grail function, 
Eq. (0), corresponds to the "naive" tunnelling factor, Eq. (P). Apparently, 



the total cross section is exponentially growing for {mw ^) \s' ^ Mq, 
but still small within the region of validity of expansion ( p^ ) . As anticipated 
above, in this energy region, it is dominated by the associated production 
of a large number of W and Z bosons. 



TT 



(nw) = 

aw 



I 6^/3 + 0(^2 



(16) 



Unfortunately, nothing is known about the behaviour of the holy grail func- 
tion for vs' around or above the barrier energy M^^^^^^^. The different terms 
in the perturbative expansion of F^ become comparable in size, and the 



perturbative expansion breaks down, just in this most interesting region. 
Unitarity and other arguments along with various assumptions have been 
used to argue [pH] that the decrease of the holy grail function may well level 
off at values of order F^^^ ~ 1/2, leading to unobservably small cross sec- 
tions of electroweak B + L violation. However, this question is not finally 
settled. 

It is very remarkable that the contribution of QCD-instantons to 
deep inelastic scattering strongly ressembles Eqs. (0)-(|lB), as first observed 
for g*g scattering in Refs. j^, ^ and elaborated for 7*51 scattering in Ref. 
O. In the Bjorken limit (c. f. Fig. for the kinematics). 
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(17) 

the total subprocess cross section for instanton-induced chirality violation is 

found to have the following structure 



where |, |, | 
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The effective instanton size /)* in Eqs. (18, T^, acting as the characteristic 
linear dimension, as well as the scaling form of F^'^^{x') are in accordance 
with the substitution rule (^ from QFD to QCD. The average gluon multi- 
plicity is found to be M, 



vr 



^ ^9 ) = — 



a. 




(20) 



It is exactly this similarity between QFD and QCD instanton-in- 
duced scattering processes which makes the study of the latter at HERA so 
interesting. 



3 Instanton-Induced Contributions to Structure 
Functions 

In this Section, let us sketch the essential steps in the pioneering calcula- 
tion of QCD-instanton contributions to the (nucleon) structure functions in 
Ref. [Q|. Along the way, we shall emphasize the basic ingredients as well 
as the inherent limitations. Finally, a state of the art evaluation of the 
/-induced contribution to the F2 structure function of the proton will be 
presented. 

First of all, it is argued [l| that the celebrated factorization theorem 
remains valid beyond conventional perturbation theory and allows to express 
the instanton contribution to the nucleon structure functions Fi^2 in the 
familiar form 

p=g,q,q -^^Bi X \ X / Y ^ J 

(21) 

In Eq. (pi|), ai = 1/2, a2(xBj) = xbj, x is the Bjorken variable of the 7*- 
parton subprocess , and fx is the factorization scale separating "hard" and 
"soft" contributions to the cross section. The distributions p{z, /x) of par- 
tons p in the nucleon absorb all information about the dynamics at large 
distances and, as usual, are to be taken from experiment. By virtue of 
Eq. (^T[), the theoretical efforts in Ref. |j|] concentrate on calculating the 

"parton" structure functions T^ (x, ^,as(At^)) in the instanton back- 
ground. For a detailed discussion on the familiar and important problem of 
infrared (IR) divergencies (associated with integrations over the instanton 
size) , we have to refer to Ref. ||l[ . In summary, it is claimed that these diver- 
gencies may be consistently absorbed into the parton distributions p{z, //), 
and an unambiguous, IR-protected contribution from small instantons may 
be isolated. 

According to the optical theorem the parton structure functions 
J-^ ^ are related to the imaginary part of the forward virtual photon-parton 
inatrix element (c.f. Fig. |2|) 

TPr°" = i/d'^e^''-(parton(p), A|r{j^(z)i,(0)}|parton(p), A). (22) 

The calculation of the instanton-induced contribution to the parton 
structure functions then involves the following steps: 



• The path integral expression for the matrix element (^) in Euclidean 
space is expanded about the instanton/anti-instanton pair configura- 
tion, defined via the so-called valley method [Q, ^P|. 

• Next, the integrations over the large number of collective coordinates 
associated with the // configuration have to be performed. 

• After Fourier transformation, the last step consists in rotating the 
result to Minkowski space and, thereafter, taking the imaginary part. 

After a long and tedious calculation, heavily exploiting the light-cone 
approximation, Balitsky and Braun |1[ succeeded in performing these steps. 
Their final answer for the instanton-induced contribution to the gluon and 
quark structure functions, derived in the Bjorken limit (c. f. Fig. g). 
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(24) 
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+ 26 



] S{C*) 



(25) 



where eg are the electric charges of the quarks, 5 = 11 — (2/3)nj, and 
d ~ 0.00363 (for n/ = 3 massless flavours) is a constant which enters the 
expression for the instanton density |^. 

The classical action of the instanton/anti-instanton pair, S{S,), is the 
most important ingredient in Eqs. (p^, [25|) , since it enters in the exponent. 

■^ For any fixed values of the collective coordinates {r}, the pair configuration (77 
"valley") is required to minimize the action within the subspace orthogonal to dG /dri. 
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a) 



b) 



Figure 2: The contribution of an // pair to to the structure function of a gluon (a) 
and of a quark (b). Sohd Unes are quark zero modes in the case that they are ending 
at the instanton (anti-instanton) , and quark propagators in the // background 
otherwise. Note that the black blobs denoted by /(/) are often referred to in the 
text as the /(/)-subprocess. 



Due to conformal invariance, it only depends on the following combination 
|0, pll of collective coordinates, 



C 



Plpj 



(26) 



with R^ = {x[ — xjf' being the instanton - anti-instanton separation, and 
piiPj their sizes, respectively. For large ^, the action ressembles a "dipole" 
form [H, I P 

(27) 



s{i) = i-^^+o{Hi)/e). 



Finally, in this "dipole" approximation, the effective conformal parameter 



11 



^=K and instanton size p* entering Eqs. (p3, 25), read 



e* 



„ l + x 

\ — X 

An 12 



(28) 
(29) 



At this point a number of important remarks should be made. 

Despite the complications associated with the 7*-parton dynam- 
ics (c.f. Fig. g), the gluon and quark structure functions ( p4|[25| ) appar- 
ently exhibit the typical signatures of an individual /-subprocess cross sec- 
tion (18, 1£). Of course, the /-subprocess variables x',Q' appearing in 
Eqs. (15, 1£) are integrated over here (c.f. also Fig. Q) and effectively sub- 



stituted by the appropriate 7*-parton variables x, Q. Let us note, in par- 
ticular, that the approximate expressions (28, 2^ agree with the solutions 
of the "saddle-point equations" - associated with the integrations over the 
collective coordinates - in case of an individual /-subprocess y, g]. 

The applicability of Eqs. (24, 25) is restricted to sufficiently large x 
(c.f. Eq. (p8|)), since their derivation was based on the large ^ ("dipole") 
approximation (p7|) for the action. A further technical requirement is 1—x ^ 



Jas{p* ), excluding the neighbourhood of x = 1. 

In general, the // interaction, Uint{(,*) = S{S,*) — 1, describes the 
emission and absorption of gluons from the instanton to the anti-instanton 
and vice versa (wavy lines between instanton and anti-instanton in Fig. ^). 
It generates via the Cutkovsky rules all final state tree-graph corrections to 
the leading semi-classical result (for a formal proof, see Ref. [^3|). These 
final state corrections are well known to exponentiate ||l^. However, it has 
been argued that some initial state and initial state - final state corrections 
exponentiate as well [EJ] and might give rise to additional corrections of 



order 47r/asC'((l 



in the exponent. 



The pre-exponential factor in Eqs. (|2j,|25D is calculated only to lead- 
ing accuracy in the strong coupling and up to corrections of order 0(1 — x). 
This is largely due to the fact, that I - I interactions have been essentially 
neglected in the prefactor, unlike the // action in the exponent. 

• On the one hand, this refers to the treatment of the "current quark" 
propagating in the // background (c.f. Fig. ^). Its presence gives 
rise to great technical complications and, correspondingly, the results 



(24, 25) only account for the first nontrivial terms in the cluster expan- 
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Figure 3: The effective conformal parameter ^* — 2 (solid) and inverse instanton 
size p^ (dashed), obtained as solutions of the saddle point equations, Eq. (pO|), for 
a range of Q values (in GeV) and nj = 3. 



sion 1 25] of the full current quark propagator in terms of the known 
propagator |26| in the background of a single (anti-)instanton. 



On the other hand, this refers to the evaluation of the functional deter- 
minants entering the pre-exponential factor. An improvement based 



on the // valley is under way [27|. 



Despite these considerable formal restrictions in the derivation of 
the gluon and quark structure functions (^, 25), it is tempting to try and 
evaluate these results within an experimentally accessible regime of xbj and 



Q . In order to hopefully enlarge the kinematical region where Eqs. (24, 25) 



may be qualitatively trusted, we have heuristically applied the following 
"improvement" steps [^ (see also Ref. jl|): 



Throughout, in Eqs. (24, 2^, we use the // valley action S'(^) of 
Refs. m, ^, rather than its "dipole" approximation (|27|). This ac- 



tion represents an extension of Eq. (27) to arbitray values of ^ > 2. It 
is worth pointing out that it continuously interpolates between S = 1 
(the sum of individual instanton and anti-instanton actions) at large ^, 
and 5 = for ^ = 2. These limiting situations correspond to a widely 
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separated, non-interacting //-pair for ^ ^ oo and the perturbative 
configuration of a collapsing and annihilating //-pair for .^ — > 2 (i.e. 
for R^ 0,pi = pj). 

We replace the approximate expressions (2^, 29) for the effective con- 
formal parameter ^* and instanton size p*, respectively, by the solu- 
tions of the exact "saddle-point" equations ||3|, ||, ||, 

Svr 



Vip* = , -u Vl7^2g^(C,), (30) 

Ois{p* ) 

Qp* = — ^(e*-2)5'(c.)-/.*^ — ^ 5(e*), 

as{p*) dp^\as{p*)J 

where 5"(^) is the derivative of the valley action 5(^) with respect to ^. 
We have numerically solved Eqs. ( |30[ ) (see Fig. ^), using the two- loop 
expression for the running coupling as(p~^) with three active flavors, 

,(3) 



and the value J^jjL = 365 MeV. It corresponds to asim-j-) = 0.33 pq| . 

The /-contribution to the nucleon structure function was finally ob- 
tained from convoluting the "improved" Eqs. (P^,|25|) with very simple phe- 
nomenological expressions ||2^ for the gluon, u quark and d quark distribu- 
tions, g{z) = (3/z)(l - zf, u{z) = (2/yi)(l - z)\ d{z) = {l/^z){l - zf, 
respectively. In view of the qualitative nature of this study, they turn out to 
be quite adequate for a factorization scale p ~ p~^, which is natural in this 
context. It turns out that over the whole xb] range considered the j*g con- 
tribution to /<2 dominates. The sea-quark contributions can be neglected 
throughout the xbj range considered. 

The resulting instanton-induced contribution to the structure func- 
tion F2 of the proton is displayed in Fig. ^. The expected very strong rise 
of the /-induced contribution with decreasing xbj is both apparent and sug- 
gestive! 

Unfortunately, any further conclusions directly reflect the (xbj,^^) 
region where the above approximations are supposed to hold. For instance, 
the dashed lines in Fig. |^ define the boundaries of various "fiducial" regions 
corresponding to values of ^(.^^(xbj, Q^)) > 0.5, 0.4, 0.3. As mentioned 
before, some authors ||2^ have advocated /Vin ~ 5'(C*)min = 1/2 as a 
saturating value for the holy grail function /^'QCD^ rpj^^ minimum value of 
Q considered in Fig. ^ is determined by the requirement that the effective 
instanton size should be sufficiently small. At Q = 10 GeV one finds /9,f ~ 1 
GeV-i (c.f. Fig. I). 
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iogio[i^i'^(^Bj,g' 




-lOoT 



Figure 4: The logarithm of the instanton-induced contribution to the structure 
function F2 of the proton, log;^Q[i^2 (^^BjiQ )], Eq. (|2l]), as a function of Xbj, 
for a range of Q values (in GeV) and Uf = 3. The curves denoted by "data" are 
to roughly represent the trend of the experimental data for F2 within the same 
range of Q. The dashed curves correspond to constant values of the valley action, 

S{UxBy,Q^))- 



4 Phenomenology of Instanton Induced Particle 
Production at HERA 

There are three main reasons which favour experimental searches for in- 
stanton-induced "footprints" in the multi-particle final state over searches 
via the structure functions, being the most inclusive observables in deep 
inelastic scattering. 

• On the one hand, the only experimental signal for QCD-instantons in 
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the structure functions could be in form of an excess over the expected 
inclusive leptoproduction rate. However, enhancements at small xbj 
are also expected from other competing mechanisms like "(perturba- 
tive) Reggeization" . Therefore, the structure functions are only of 
limited value in searches for manifestations of QCD-instantons. 

• On the other hand, as we shall see, the instanton-induced final state 
is distinguished by a quite spectacular event topology together with a 
characteristic flow of flavour quantum numbers. 

• Furthermore, the additional possibility of imposing experimental cuts 
on kinematical variables of the final state may well allow to restrict 
the I-subprocess variables x' ,Q', . . . within a theoretically controllable 
regime despite small xbj- Along these lines one may hope to bridge 
the substantial gap between the regime of larger x' >O{0.1), where 
the I-subprocess cross sections may be theoretically estimated, and 
the small xbj regime, xbj <O(10~^), where the bulk of HERA data is 
accumulating at present. 

To elaborate on the last two aspects is the purpose of this Section. 

A graphical display of the modulus squared of the relevant j*g ma- 
trix element, along with various four- momenta of interest, is presented in 
Fig. ^. Its structure suggests that the instanton-induced contribution to the 
differential cross sections for j*g{q) — > • • • can be written in a "canonical" 
convolution form, as familiar from perturbative QCD: 

(31) 

The integrations in Eq. ( [3l| ) extend over the variables x' , Q'^, referring as 
in Sect. 2 to the instanton-induced subprocess (denoted by / in Fig. ||), 

s' = {q'+pf 2pq' s' + Q'^- ^ ^ 

Their definition is completely analogous to the standard e^A^ variables Q^ = 
—q^ and xbj = Q^ /2Pq referring to the nucleon target of momentum P. The 
parton (gluon) momentum fraction with respect to the proton is z = xbj/x = 
{pq)/{Pq) and 

< XBj < X < x' < 1. (33) 
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Figure 5: Graphical display of the instanton-induced contribution to the cross 
section of ^*g scattering for AQ = 2nj = 6. 



The conditions for the vahdity of Eq. ( |3l| ) (beyond the set of approx- 
imations inherent in Ref. ll|), as well as a determination of the "splitting 
function" f^J ,{z', Q'"^), associated with the propagation of the current quark 



in the instanton background, are presently under active investigation |3C]. 
To establish a structure of type ( |3l] ) is quite an important task both from 
a theoretical point of view and also for further studies of instanton-induced 



phenomenology by means of Monte Carlo methods [31|. For the time being, 
we shall simply assume that Eq. (^l|) is valid approximately. 
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It is quite plausible that the "splitting function" f!^J,{z',Q''^) in 
Eq. ( [31] ) only exhibits a relatively mild dependence on its variables. In 
contrast, the I-subprocess cross sections da , ix' , Q'^, . . .) bring in the main 
dependences (of exponential type) and, of course, are most interesting from 
the physics point of view. Accordingly, we have concentrated our theoretical 
efforts in Ref. ^ on calculating the crucial observables characterizing the 
I-subprocess, such as normalized 1,2,.. .-parton inclusive cross sections, 

;^r^<^'^' «'■■■■'• '^^' 

along with the respective average parton multiplicities, transverse momen- 
tum (flow), etc. Corresponding to the restrictions discussed in Sect. 3 in the 
context of the structure functions, the I-subprocess variables should not be 
too small, x'>C'(0.1), Q'>C'(10 GeV), say. 

Since a calculation of the "splitting function" in Eq. (|3l| ) is still under 
way |3C], a discussion of expected event rates has to be deferred to a later 



stage. In the present analysis, we only make use of information abstracted 
from our calculations of the I-subprocess observables Q along with HERA 
kinematics. This is sufficient, however, to obtain important insight into 
the expected event topology in the most interesting regime of small xbj- 
Moreover, the connection between kinematical quantities measurable in the 
laboratory system and the variables controlling the instanton subprocess 
may be studied. 

Specifically, we use the following set of working hypotheses about 
the I-subprocess: 

i) Isotropy: In its cm. system, q'+p = 0, the instanton-induced multi-parton 
production is supposed to proceed isotropically. We may imagine a "fireball" 
in S-wave configuration, decaying into gluons and at least 2nf — 1 quarks, 
including strangeness (!) and possibly charm, if kinematically allowed (c. f. 
Fig. I). 

ii) Dependence on x', Q'^: The I-subprocess cross sections da , (x' , Q'^ , . . .) 
are expected to strongly decrease with increasing Q'^ for fixed x' and to 
strongly increase with decreasing x' for fixed Q'^ (c. f. Sect. 2). As discussed 
above and in Sect. 3, it remains uncertain, however, how long the cross 
sections continue to increase towards x' -^ 0. 
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Hi) Multiplicity: The total multiplicity associated with the I-subprocess is 
expected to be large, 

(ng+q(x',g'2)) ^ o (^^) + 2nf - 1 >O(10), (35) 

on the parton level, leading typically to 0(20-^30) particles after hadroniza- 
tion. In Fig. |^, we display {ug^qix' , Q'^)) as calculated in Ref. Q. 

At small values of x' and large Q'^, the multiplicity obtains a large 
contribution from gluons and peaks around x' ~ 0.2 H- 0.3, whereas at large 
x' the 2nf — 1 = 5 produced quarks dominate. The peaking of the gluon 
multiplicity at small, non-vanishing x' actually has an appealing interpreta- 
tion: 

For large values of Q'^ the coefficient of 7r/2as in Eq. p5| ) turns 
out to be [P 4(^^, — 2)S"(^*), involving the derivative of the valley ac- 
tion with respect to the conformal parameter S,, taken at the saddle point 
value S.^ix^Q''^) (c.f. Fig. ^. As discussed in Sect. 3, the full valley ac- 
tion smoothly interpolates between a non-interacting, infinitely separated 
instanton/anti-instanton pair for ^ — > oo (probed for x' ^ 1) and the per- 
turbative vacuum for .^ ^ 2 (probed for x' -^ 0). Hence, in both limits 
a decrease of the gluon multiplicity matches well the intuition! The peak 
of the multiplicity inbetween corresponds to the maximal variation of the 
action with ^. 

We also note the substantial increase of the gluon multiplicity with 
increasing Q'^, which at large Q'^ mainly reflects the running of a^ in Eq. 

(D. 

iv) k± signature: The transverse momenta of the partons emerging from the 
instanton subprocess and the one of the current-quark jet (c.f. Fig. |5|) are 
expected to be "semi hard" , typically of order |^ |l| 

I, I < ^±tot > IT Vs' Vs' . ^ J , .o«^ 

\k\i\ ~ = ~ Oo ; I C 1, and (6b j 

'' <n> 4<n> 2 ^ ' 



I A;± current quark I ~ yOigQ'^. (37) 

Given this plausible generic input i) - iv), we may now ask, how instanton- 
induced events would look like in the Hl/ZEUS detectors. 

First of all, we observe that as a direct consequence of the isotropy 
assumption i), the (pseudo) rapidity distribution of a single final state parton 

19 



[n,+,{x',Q''')) 




Figure 6: The average, total parton multiplicity associated with the /-subprocess 
as a function of x' , for different values of Q' and n/ = 3, from Ref. Q. 



in the I-c.m. system takes the form 

1 dcT(^)(x',Q'2,r/7) _ 1 <n(x',Q'2)> 



a 



{I ) tot 



dr]id(j)i 



(38) 



i.e. it is strongly peaked in (pseudo) rapidity rji 
r]j = with a half width of 

Ar]i w ±0.9. 



47r cosh (r//) 2 

— lntan(0//2) around 



(39) 



The shape and width of the distribution in pseudo rapidity remains, of 
course, very similar in the HERA laboratory system, for kinematical con- 
figurations where the /-cm. system is dominantly boosted longitudinally. 
Depending on the values of the various subprocess variables, the peak posi- 
tion {r]j = 0) fluctuates in general over the available range of rj (for given 
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rEBj and y viz. Q"^): 



.lab 



ryft^ = -In 



1 1. ( ^^ ^Bj 



yEe X^ 



+ x^(l - y) + (1 - x^)(- - 1) 

X X 



-2y (i-2/)2;7(i-x-y)(— -i)cos X 



(40) 



In analogy to the standard y variable, we have introduced in Eq. (|4^) the q' 

momentum fraction 

o'v 
< x^ = -i-^ < 1. (41) 

qp 

The variable x denotes the azimuthal angle of the vector q' in the HERA lab- 
oratory frame. Due to momentum conservation, \q[_\ sin^ = —^i/ current quark) 
the component of the current-quark momentum out of the ee'P scattering 
plane. 

In Fig. 1^, the dominant matrix element and {rj , i?i> )-plot of a typ- 
ical instanton-induced event is contrasted with the canonical two-jet config- 
uration in perturbative QCD. 

Clearly, the "0th level" signature to watch out for is a densely pop- 
ulated hadronic "band" in the (r/^^^, c/)^'^^) -plane, centered at a fluctuating 
value (|4^) of rjj . This striking multi-hadron final state originates from 
8 H- 10 "semi-hard" jets (c. f. Fig. |6|), always includes strangeness and is 
characterized by a width Arjj = ±0.9. Let us point out two observables, 
which - on an event-by-event basis - appear to be particularly sensitive to 
this event structure. 

• The (transverse) energy flow, dEn\/dr] (integrated over (j) ), will 
exhibit a strong enhancement at the position r] = rjf of the "band" , 
since each of the 8-^10 instanton-induced jets contributes a comparable 
energy into a single rj bin of width ~ 1.8. If, in addition, the current- 
quark jet is isolated from the "band" (see below), one even expects a 
double-peak structure in dEn\/dr] . The energy flow signature may 
well be less affected by hadronization than patterns associated with 
individual tracks. 

• Pseudo sphericity: The usual event-shape variables like sphericity and 
aplanarity should be useful tools in analysing the manifestations of 
an isotropic instanton-induced subprocess in the final state. Of par- 
ticular sensitivity appears the so-called pseudo sphericity [^], which 
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Figure 7: The dominant matrix element and [t] , (jr) -plot of a typical 
instanton-induccd event is contrasted with the canonical two-jet configuration in 
perturbative QCD. 
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incorporates only transverse information from the event in terms of 
the azimuthal angles (pi of the N final state hadrons: 



Pseudo sphericity = 1 — — : 



\ 



N N 



(^cos(/),)2 + (^sin,/),)2. (42) 



1=1 i=l 



Apparently, it equals to 1 in the cm. system of an isotropical event 
and vanishes for a single "collimated" jet. 

As a next level of sophistication, we study the effects of kinemati- 
cal cuts on suitable final-state variables with the aim, to further enhance 
the event topology and to unfold and/or restrict the /-subprocess variables 
x' , Q'. Of course, this has to be achieved 

• without affecting significantly the expected size of the /-subprocess 
cross section do"'-^ (by exploring input ii) above); 

• such that the /-subprocess is invoked in a kinematical region of {x' , Q'), 
where it induces a high average multiplicity (ng+q) according to Fig. ^. 

The key aspect is to focus on an event topology corresponding to an 
isolated (semi-hard) current-quark jet (c. f. input iv) above: k± ^ \/(xJ^) 
in addition to a hadronic "band" in [t]^^^ , (p^'^^) as discussed in "level 0" 
above. To this end, let us consider the constraints on the (internal) subpro- 
cess variables x^, x' , Q' . . ., implied by the following reasonable separation 
criteria: 

• The hadronic "band" should be contained in the (central) detector, 
such that its peak position satisfies |??/''^|<1, say (in practice, this 
upper bound for iiiY^ may well be pushed up to ~ 2). 

• The current-quark jet is required to be separated in r/ from the 
hadronic "band" (centered around 77^*^^) by 

An = rv'^l^ - r)!*^^ i ~ ~ ^'^ °^ (A-i] 

^7 — 7/ '/current quark S > _j_ 1 5 ^ ' 

• A minimal transverse momentum, A;_l current quark ^4 GeV, is required 
for the current-quark jet. 
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In Fig. ^ we have displayed the resulting restrictions on the internal 
subprocess variables x^, x' , for a typical set of fixed "external" parame- 
ters XBj,y,a^- Apparently, after imposing the isolation requirements for the 
current-quark jet we are left with two allowed, 'triangular' regions 1 and 2, 
in the (3;^,x')-plane. The central portion of Fig. ^ is excluded by the cut 
on At/ (solid lines), with the left- (right-) hand boundary corresponding to 
At/ = +1.5 (—1.5). The excluded portion on the left (short dashes) refers 
to the hadronic "band" being centered within 1 < "rff^ < 2.2, with rjY^ = 1 
located on the right. Hence, if a value of ifj above 1 is experimentally tol- 
erable, the allowed region 1 increases significantly. Finally, the main effect 
of the requirement k± current quark > 4 GeV , is to set a lower limit to the 
involved values of x' , and to exclude the region x^ very close to 1 (long 
dashes) . 

These results demonstrate that, indeed, kinematical cuts of the type 
considered here, may well restrict the /-subprocess variables x' , Q' to regions 
where the computation Q of dcrL (x', Q') may be trusted (e.g. within regime 
2 of Fig. I). 

According to our input ii) above and Fig. P, regime 1 in Fig. |8| will 
presumably be associated with considerably higher rates, since it typically 
corresponds to much smaller values of Q' and x' than regime 2. From Fig. y 
we also infer a comfortably high total parton multiplicity (Tig+q) ~ 8 in 
regime 1, as well as an energy/parton = \/s'/8>3 GeV in the /-cm. system. 

Let us finally illustrate in Figs. |l^, |l^ an event in the HERA labo- 
ratory frame, corresponding both to a striking signature (isolated current- 
quark jet along with a densely populated hadronic "band") and favourable 
rate/multiplicity conditions for the instanton subprocess. We note the corre- 
sponding event-shape variables (on the parton level in the laboratory frame) 
as calculated from the instanton-induced partons within the "band" 

Sphericity|[g^ ^ 0.45, 

Aplanarityfi^jj, « 0.22, (44) 

Pseudo sphericityLj, ~ 0.90. 

Of course, due to the isotropy of the instanton-induced subprocess, the 
sphericity, aplanarity and pseudo-sphericity variables essentially adopt their 
maximal values of 1, 1/2 and 1, respectively, in the /-cm. frame. Ap- 
parently, the pseudo sphericity is least affected by the Lorentz transfor- 
mation into the laboratory frame and continues to reflect the underlying 
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Figure 8: In the displayed regions 1,2 of the (a;-y,x')-plane, the current-quark 
jet is isolated from the instanton-induced "band" (centered at ijf" ). The excluded 
domains are: -1.5 < Ar] = Tfff° - rf^^„^^^ qu^rk ^ 1-5 (sohd hues), 1 < rff" 
(short dashes) and A;_|_ current quark ^ 4 GeV (long dashes). The parameters are 
Ep = 820 GeV, Ee = 30 GeV, xbj = 0.001, y = 0.5 (Q = 7 GeV), x = 0.025 
(Vs = 44 GeV) and q' azimuthal angle X = 0- 



isotropy. Clearly, for a quantitative discussion, hadronization effects have 
to be included |3^, which may well wash out somewhat the traces of the 
underlying isotropy in the various event shape parameters. 



5 Summary and Outlook 

The search for QCD-instanton induced events at HERA is well worth the 
effort: 

First of all, these "anomalous" processes are predicted to occur 
within standard QCD. Secondly, there is a close analogy to electroweak B+L 



25 




Figure 9: Lines of constant invariant mass of the I-subprocess, v s', and constant 
Q' (dashed), versus the subprocess variables x' and x^, for the same parameters as 
in Fig. y. The numbers at the crossing points denote the respective /-subprocess 
multiphcities (?^g+q) according to Fig. ^. 



violating processes, as was discussed in detail in Sect. 2. While a promising 
search for anomalous electroweak events is only possible in the far future, 
the search for manifestations of QCD-instantons can start right now, in deep 
inelastic e^p scattering at HERA. 

Besides summarizing the essence and limitations of the theoretical 
calculations involved |j|], we have presented in Sect. 3 a state of the art evalu- 
ation of the instanton-induced contribution to the nucleon structure function 
i^2- It rises strongly with decreasing xbj and tends to reach the size of the 
experimental data around xbj ~ 0.1 -^ 0.25. Unfortunately, due to inher- 
ent uncertainties, the calculation cannot be trusted anymore for xbj <0.35, 
say. Nevertheless, the trend is very suggestive! However, enhancements of 
the inclusive leptoproduction rate at small xbj are also expected from other 
competing mechanisms like "(perturbative) Reggeization" . Therefore, the 
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Figure 10: 3d-niomcntum display (in GeV) for a typical instanton-induced event 
before hadronization, in the HERA laboratory frame, satisfying the kinematical 
cuts discussed in the text. The current-quark jet is well isolated. Not shown are 
the incoming proton {-\-z direction) and e {—z direction), as well as the proton 
fragments. The parameters are as in Fig. g and moreover, x^ = 0.32, x' = 0.1, 
such that Q' ^8 GeV, \/s' « 24 GeV and (ng+q) ^ 8. 



structure functions appear only of limited value in searches for "footprints" 
of QCD-instantons. 

In Section 4 we have reported first phenomenological results of our 
ongoing broad and systematic investigation of the QCD-instanton induced 
hadronic final state Q. Since a calculation of the "splitting function" - 
associated with the propagation of the current quark in the instanton back- 
ground - is still in progress, a discussion of expected event rates has to 
be deferred to a later stage. In the present analysis we only made use of 
information abstracted from our calculation of the instanton-subprocess ob- 
servables 0] along with HERA kinematics. This was sufficient, however, 
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Figure 11: Lego plot of the transverse energy in GeV (before hadronization) 
for the same instanton-induced event as in Fig. nfl, satisfying the kinematical cuts 
discussed in the text. Not shown are the scattered e and the proton fragments. 



to obtain important insight into the expected event topology in the most 
interesting regime of small xbj • 

• The "0th level" signature to watch out for is a densely populated 
hadronic "band" in the [r]^^^, (j)^''^^) -plane (c. f. Fig. ^), centered at some 
fluctuating value r]j . This striking multi-hadron final state originates 
from 8-^ 10 "semi-hard" jets (c. f. Fig. |6|), always includes strangeness 
and is characterized by a width Arjj = ±0.9. It directly reflects 
an underlying instanton-subprocess, associated with the formation of 
an S'-wave "fireball", which then decays isotropically into gluons and 
at least 2nf — 1 quarks. Observables which are particularly sensitive 
to this event structure are e.g. the (transverse) energy flow and the 
so-called pseudo sphericity. The energy flow signature may well be 
less affected by hadronization than patterns associated with individual 
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tracks. 

• As a next level of sophistication, we have studied kinematical cuts on 
suitable final-state variables, which help to further enhance the event 
topolology and to unfold and/or restrict the (Bjorken) variables of the 
instanton subprocess within a theoretically controllable regime, despite 
small XBj- Along these lines one may hope to bridge the substantial 
gap between the kinematical region, where the instanton-subprocess 
cross sections may be theoretically estimated [Q], and small xbj val- 
ues, XBj <O(10^'^), where the bulk of HERA data is accumulating at 
present. 

In summary, experimental searches for instanton "footprints" in the multi- 
particle final state appear to be much more promising than searches via the 
structure functions. 

Finally, let us briefly mention some related theoretical and phe- 
nomenological issues presently under study. Theoretical work is in progress 
to improve the pre-exponential factors, affecting quite strongly the predic- 
tions for structure functions and the various subprocess cross sections. This 
refers in particular to a more reliable evaluation of the functional deter- 
minants [^] in the instanton/anti-instanton valley background beyond the 
dilute instanton-gas approximation. 

Of great importance for further studies of QCD-instanton phenome- 
nology is the task of establishing a convolution form [^] of the 7*-parton 
multi-particle cross sections in terms of "splitting functions" and instanton- 
subprocess cross sections (c. f. Eq. (^)). Once the "splitting functions" 
have been isolated and calculated, we hope to come forward with predic- 
tions for the rate of instanton-induced multi-particle events. We are then 
ready to study the instanton-induced multi-particle final state by means of 
a Monte Carlo based event generator p^]. Only after including effects of 
hadronization and background will it be possible to address the crucial ques- 
tion: How many "anomalous" events are needed to establish the "discovery" 
of an instanton at HERA? 
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